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1. INTRODUCTION 



These four lectures are devoted to an elementary introduction to Quantum Chromo-Dynamics (QCD), 
the theory of strong interactions. Four lectures are not much. So after a general introduction I will 



concentrate on the basic principles and the main applications of perturbative QCD (for reviews of the 
subject see for example |JJ|, [Q], [J3j] ) - I will try to put the main emphasis on ideas with only a minimum 
of technicalities. 



<D 

At present most of particle physics is well understood in terms of the Standard Model (SM), which 
is a gauge theory of strong and electroweak interactions based on the group SU (3) SU(2) U (l)y . 
The 577(3) factor is the colour group of QCD, while 5*7(2) U(l)y is the Glashow-Weinberg-Salam 
electroweak symmetry group. The electroweak symmetry is spontaneously broken down to U{\)q, the 
phase group of the electric charge Q e (that is, Q is the charge operator in units of e, the proton charge), 
different from the U (l)y of weak hypercharge: Q = t$ + Y/2, where ts is the third component of the 
weak isospin generator of SU (2). The group SU (3) U (1)q is believed to be an exact gauge symmetry 
of nature. The corresponding gauge bosons, the eight gluons and the photon are massless. Matter fields 
include three families of coloured quarks and colourless leptons. Quarks and gluons (q and g) are the 
fields that have strong interactions described by QCD. The statement that QCD is a gauge theory based 
on the group 5*7(3) with colour triplet quark matter fields fixes the QCD lagrangian density to be (a 
summary of the general formalism of gauge theories is presented in Appendix): 

i 8 n f 

A=l j=l 

Here: qj are the quark fields (of nf different flavours) with mass mf, I) = D^j 11 , where 7^ are the Dirac 
matrices and is the co variant derivative: 

13/1 = 9^- iesgf,; (2) 

e s is the gauge coupling, later we will mostly use, in analogy with QED 



g M = J^a t 9a wh ere -4 = 1,8, are the gluon fields and t A are the SU(3) group generators in the 
triplet representation of quarks (i.e. tA are 3x3 matrices acting on q); the generators obey the commuta- 
tion relations [t , t B ] = iCABct C where Cabc are the complete antisymmetric structure constants of 
SU{3) (the normalisation of C A bc and of e s is specified by Tr[t A t B ] = 5 AB /2); 

< = d a g A -d v g A - e s C A BC9*9 C v (4) 

The Feynman rules of QCD are listed in Fig 1 . The physical vertices in QCD include the gluon- 
quark-antiquark vertex, analogous to the QED photon-fermion-antifermion coupling, but also the 3-gluon 
and 4-gluon vertices, of order e s and e 2 respectively, which have no analogue in an abelian theory like 
QED. In fact the QED lagrangian density is given by: 

C = -\f^F^ + ^(i^-m^ (5) 

with the covariant derivative given in terms of the photon field 4 M and the charge operator Q by: 

D u = d u - ieA^Q (6) 

and 

iV = 9„A U - d v A^ (7) 

In QED the photon is coupled to all electrically charged particles but itself is neutral. In QCD the gluons 
are coloured hence self-coupled. This is reflected in the fact that in QED F^ is linear in the gauge field, 
so that the term F 2 in the lagrangian is a pure kinetic term, while in QCD F A V is quadratic in the gauge 
field so that in F A2 we find cubic and quartic vertices beyond the kinetic term. Also instructive is to 
consider the case of scalar QED: 

C = -\F^F UV + ( J D At 0)t( J D^)-m 2 (0t0) (8 ) 

For Q = 1 we have: 

{D^(D^) = [d^)\d^) + ieA^[(d^)U ~ <t> j (d^)} + e 2 A^<pU (9) 

We see that for a charged boson in QED, given that the kinetic term for bosons is quadratic in the 
derivative, there is a two-gauge vertex of order e 2 . Thus in QCD the 3-gluon vertex is there because the 
gluon is coloured and the 4-gluon vertex because the gluon is a boson. 

The QCD lagrangian in eq.(|]) has a simple structure but a very rich dynamical content. It gives 
rise to a complex spectrum of hadrons, it implies the striking properties of confinement and asymp- 
totic freedom, is endowed with an approximate chiral symmetry which is spontaneously broken, has a 
highly non trivial topological vacuum structure (instantons, U(1)a symmetry breaking, strong CP vio- 
lation (?)...), an intriguing phase transition diagram (colour deconfinement, quark-gluon plasma, chiral 
symmetry restauration, colour superconductivity, ...). 

Confinement is the property that no isolated coloured charge can exist but only colour singlet 
particles. For example, the potential between a quark and an antiquark has been studied on the lattice 
and it has a Coulomb part at short distances and a linearly rising term at long distances: 

Vm „ C F p& +.... + or] (10) 
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Fig. 1: Feynman rules for QCD. The solid lines represent the femions, the curly lines the gluons, and the dotted lines represent 
the ghosts. 



with Nc the number of colours (Nc = 3 in QCD). The scale dependence of a s (the distance r is Fourier- 
conjugated to momentum transfer) will be explained in detail in this course. The linearly rising term 
makes it energetically impossible to separate a q — q pair. If the pair is created at one space-time point, 
for example in e + e~ annihilation, and then the quark and the antiquark start moving away from each 
other in the center of mass frame, it soon becomes energetically favourable to create additional pairs, 
smoothly distributed in rapidity between the two leading charges, which neutralise colour and allow the 
final state to be reorganised into two jets of colourless hadrons, that communicate in the central region by 
a number of "wee" hadrons with small energy. It is just like the familiar example of the broken magnet: 
if you try to isolate a magnetic pole by stretching a dipole, the magnet breaks down and two new poles 
appear at the breaking point. Very often in QCD one computes inclusive rates for partons (the fields in 
the lagrangian, that is, in QCD, quarks and gluons) and takes them as equal to rates for hadrons. Partons 
and hadrons are considered as two equivalent sets of complete states. This is called "global duality" and 
it is rather safe in the rare instance of a totally inclusive final state. It is less so for distributions, like 
distributions in the invariant mass M ("local duality") where it can be reliable only if smeared over a 
sufficiently wide bin in M. 

Confinement is essential to explain why nuclear forces have very short range while massless gluon 
exchange would be long range. Nucleons are colour singlets and they cannot exchange colour octet 
gluons but only colourless states. The lightest colour singlet hadronic particles are pions. So the range 
of nuclear forces is fixed by the pion mass r ~ m^ 1 ~ 10 -13 cm : V « exp(— m n r)/r. 

Why SU(Nc = 3) co i our l The selection of SU(3) as colour gauge group is unique in view of 
a number of constraints, (a) The group must admit complex representations because it must be able to 
distinguish a quark from an antiquark. In fact there are meson states made up of qq but not analogous qq 
bound states. Among simple groups this restricts the choice to SU(N) with N > 3, SO(AN + 2) with 
N > 2 (taking into account that SO(6) has the same algebra as SU(A)) and E(6). (b) The group must 
admit a completely antisymmetric colour singlet baryon made up of 3 quarks: qqq. From the study of 
hadron spectroscopy we know that the low lying baryons, completing an octet and a decuplet of (flavour) 
SU(3) (the approximate symmetry that rotate the 3 light quarks u, d and s), are made up of three quarks 
and are colour singlets. The qqq wave function must be completely antisymmetric in colour in order to 
agree with Fermi statistics. Indeed if we consider, for example, a N* ++ with spin z-component +3/2, 
this is made up of (u "ft u "fl* u fl*) in an s-state. Thus its wave function is totally symmetric in space, 
spin and flavour so that complete antisymmetry in colour is required by Fermi statistics. In QCD this 
requirement is very simply satisfied by e a b c q a q b q c where a, b, c are SU(3) co i our indices, (c) The choice 
of SU(Nc = 3) co i our is confirmed by many processes that directly measure Nc- Some examples are 
listed here. The total rate for hadronic production in e + e _ annihilation is linear in Nc- Precisely if we 
consider R = a{e + e~~ — > hadrons) /a po i n t(e + e~ — ► fi + ' n~) above bb threshold and below mz and we 
neglect small computable radiative corrections (that will be discussed later) we have a sum of individual 
contributions (proportional to Q 2 ) from qq final states with q = u, c, d, s, b: 

R ~ N c [2.\+3.\] ~N C 1 ± (12) 

The data neatly indicate Nc = 3 as seen from Fig. 2 [HJ. The slight excess of the data with respect to 
the value 1 1/3 is due to the QCD radiative corrections. Similarly we can consider the branching ratio 
B(W~ — > e~u), again in Born approximation. The possible fermion-antifermion (//) final states are 
for / = e~, T~,d,s (there is no / = b because the top quark is too heavy for bt to occur). Each 
channel gives the same contribution, except that for quarks we have Nq colours: 

B(w ~- " irnsvc <13) 

For Nc = 3 we obtain B = 11% and the experimental number is B = 10.7%. Another analogous 
example is the branching ratio B(t~ — > e~v e v T ). From the final state channels with / = e~, u~ , d we 
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find 

S(r- - e"^) « (14) 

For Nq = 3 we obtain B = 20% and the experimental number is B = 18% (the less accuracy in this 
case is explained by the larger radiative and phase-space corrections because the mass of r~ is much 
smaller than my/)- An important process that is quadratic in Nc is the rate r(-7r° — ► 27). This rate can 
be reliably calculated from a solid theorem in field theory which has to do with the chiral anomaly: 

IV - 27) « (^) 2 ^| = (7-73 ± 0.04)(^)= eV (15) 

where the prediction is obtained for f w = (130.7 ± 0.37) MeV. The experimental result is T = (7.7 ± 
0.5) eV in remarkable agreement with iVc = 3. There are many more experimental confirmations that 
Nc = 3: for example the rate for Drell-Yan processes (see section 5.4 ) is inversely proportional to Nc- 

How do we get testable predictions from QCD? On the one hand there are non perturbative meth- 
ods. The most important at present is the technique of lattice simulations: it is based on first principles, 
it has produced very valuable results on confinement, phase transitions, bound states, hadronic matrix 
elements and so on, and it is by now an established basic tool. The main limitation is computing power 
and therefore there is continuous progress and a lot of good perspectives for the future. Another class of 
approaches is based on effective lagrangians which provide simpler approximations than the full theory, 
valid in some definite domain of physical conditions. Chiral lagrangians are based on soft pion theorems 
and are valid for suitable processes at energies below IGeV. Heavy quark effective theories are obtained 
from expanding in inverse powers of the heavy quark mass and are mainly important for the study of b 
and, to less accuracy, c decays. The approach of QCD sum rules has led to interesting results but appears 
to offer not much potential for further development. Similarly specific potential models for quarkonium 
have a limited range of application. On the other hand, the perturbative approach, based on asymptotic 
freedom, still remains the main quantitative connection to experiment, due to its wide range of applica- 
bility to all sorts of "hard" processes. Perturbative QCD will be the main subject in the following. I will 
discuss its foundations and main applications in the next sections. 



2. Massless QCD and Scale Invariance 

The QCD lagrangian in eq.(|]) only specifies the theory at the classical level. The procedure for quanti- 
sation of gauge theories involves a number of complications that arise from the fact that not all degrees 
of freedom of gauge fields are physical because of the constraints from gauge invariance which can be 
used to eliminate the dependent variables. This is already true for abelian theories and we are familiar 
with the QED case. One introduces a gauge fixing term (an additional term in the lagrangian density that 
acts as a Lagrange multiplier in the action extremisation). One can choose to preserve manifest Lorentz 
invariance. In this case, one adopts a covariant gauge, like the Lorentz gauge, and in QED one proceeds 
according to the formalism of Gupta-Bleuler. Or one can give up explicit formal covariance and work in 
a non covariant gauge, like the Coulomb or the axial gauges, and only quantise the physical degrees of 
freedom (the transverse components of the photon field). While this is all for an abelian gauge theory, in 
the non-abelian case some additional complications arise, in particular the necessity to introduce ghosts 
for the formulation of Feynman rules. There are in general as many ghost fields as gauge bosons and 
they appear in the form of a transformation Jacobian in the Feynman diagram functional integral. Ghosts 
only propagate in closed loops and their vertices with gluons can be included as additional terms in the 
lagrangian density which are fixed once the gauge fixing terms and their infinitesimal gauge transforma- 
tions are specified. We skip the detailed derivation of the complete Feynman rules in a given gauge as 
they appear in Fig 1. 

Once the Feynman rules are derived we have a formal perturbative expansion but loop diagrams 
generate infinities. First a regularisation must be introduced, compatible with gauge symmetry and 
Lorentz invariance. This is possible in QCD. In principle one can introduce a cut-off A (with dimen- 
sions of energy), for example, a' la Pauli-Villars. But at present the universally adopted regularisation 
procedure is dimensional regularisation that we will briefly describe later on. After regularisation the 
next step is renormalisation. In a renormalisable theory (like for all gauge theories in 4 spacetime di- 
mensions and for QCD in particular) the dependence on the cutoff can be completely reabsorbed in a 
redefinition of particle masses, of gauge coupling(s) and of wave function normalisations. After renor- 
malisation is achieved the perturbative definition of the quantum theory that corresponds to a classical 
lagrangian like in eq.(|IJ) is completed. 

In the QCD Lagrangian of eq.([I]) quark masses are the only parameters with physical dimensions 
(we work in the natural system of units h = c = 1). Naively we would expect that massless QCD is 
scale invariant. This is actually true at the classical level. Scale invariance implies that dimensionless 
observables should not depend on the absolute scale of energy but only on ratios of energy-dimensional 
variables. The massless limit should be relevant for the asymptotic large energy limit of processes which 
are non singular for m — > 0. 

The naive expectation that massless QCD should be scale invariant is false in the quantum theory. 
The scale symmetry of the classical theory is unavoidably destroyed by the regularisation and renormal- 
isation procedure which introduce a dimensional parameter in the quantum version of the theory. When 
a symmetry of the classical theory is necessarily destroyed by quantisation, regularisation and renormal- 
isation one talks of an "anomaly". So, in this sense, scale invariance in massless QCD is anomalous. 

While massless QCD is finally not scale invariant, the departures from scaling are asymptotically 
small, logarithmic and computable. In massive QCD there are additional mass corrections suppressed by 
powers of m over the energy scale (for non singular processes in the limit m — > 0). At the parton level 
(q and g) we can conceive to apply the asymptotics from massless QCD to processes and observables 
(we use the word "processes" for both) with the following properties ("hard processes"), (a) All relevant 
energy variables must be large: 

Ei = ZiQ, Q » rrij; Zf. scaling variables o(l) (16) 

(b) There should be no infrared singularities (one talks of "infrared safe" processes), (c) The processes 
concerned must be finite for m — > (no mass singularities). To possibly satisfy these criteria processes 



must be as "inclusive" as possible: one should include all final states with massless gluon emission and 
add all mass degenerate final states (given that quarks are massless also q — q pairs can be massless if 
"collinear", that is moving together in the same direction at the common speed of light). 

Let us discuss more in detail infrared and collinear safety. Consider, for example, a quark virtual 
line that ends up into a real quark plus a real gluon (Fig. 3). 




Fig. 3: 



For the propagator we have: 

Pr ° Pagat0r = (p + k^-m^ = 2^k) = 2E^ • 1- please (17) 

Since the gluon is massless, can vanish and this corresponds to an infrared singularity. Remember 
that we have to take the square of the amplitude and integrate over the final state phase space, or all 
together, dE^jE^. So we get 1/E% from the squared amplitude and d 3 k/Ek ~ E^dE^ from the phase 
space. Also, for m — > 0, (3 p = ^Jl — m 2 / E 2 — > 1 and (1 — p cos 9) vanishes at cos 9 = 1. This leads 
to a collinear mass singularity. 

There are two very important theorems on infrared and mass singularities. The first one is the 
Bloch-Nordsieck theorem: infrared singularities cancel between real and virtual diagrams (see Fig. 4) 
when all resolution indistinguishable final states are added up. For example, for each real detector there 
is a minimum energy of gluon radiation that can be detected. For the cancellation of infrared diver- 
gences, one should add all possible gluon emission with a total energy below the detectable minimum. 
The second one is the Kinoshita-Lee-Nauenberg theorem: mass singularities connected with an external 
particle of mass m are canceled if all degenerate states (that is with the same mass) are summed up. That 
is for a final state particle of mass m we should add all final states that in the limit m — ► have the same 
mass, also including gluons and massless pairs. If a completely inclusive final state is taken, only the 
mass singularities from the initial state particles remain (we shall see that they will be absorbed inside 
the non perturbative partem densities, which are probability densities of finding the given partem in the 
initial hadron). 

Hard processes to which the massless QCD asymptotics can possibly apply must be infrared and 
collinear safe, that is must satisfy the requirements from the Bloch-Nordsieck and the Kinoshita-Lee- 
Nauenberg theorems. We give now some examples of important hard processes. One of the sim- 
plest hard processes is the totally inclusive cross section for hadron production in e + e~ annihilation, 
Fig. 5, parameterised in terms of the already mentioned dimensionless observable R = a{e + e~ — > 
hadrons) / a po int(e + e~ — > fi + fi~). The pointlike cross section in the denominator is given by a po i n t = 
4ira 2 /3s, where s = Q 2 = 4E 2 is the squared total center of mass energy and Q is the mass of the 
exchanged virtual gauge boson. At parton level the final state is (qq + n g + n' q'q') and n and n' 
are limited at each order of perturbation theory. It is assumed that the conversion of partons into hadrons 
does not affect the rate (it happens with probability 1). We have already mentioned that in order for this 
to be true within a given accuracy an averaging over a sufficiently large bin of Q must be understood. 
The binning width is larger in the vicinity of thresholds: for example when one goes across the charm 



Fig. 4: 




cc threshold the physical cross-section shows resonance bumps which are absent in the smooth partonic 
counterpart which however gives an average of the cross-section. 

A very important class of hard processes is Deep Inelastic Scattering (DIS) 

l + N^l' + X / = e ± ,/i ± ,^,P (18) 

which has played and still plays a very important role for our understanding of QCD and of nucleon 
structure. For the processes in eq.(18), shown in Fig 6, we have, in the lab system where the nucleon of 
mass m is at rest: 

O 2 

Q 2 = -q 2 = -{k-k') 2 = AEE' sin 2 6/2; mv = (p.q); x = (19) 

2mu 

In this case the virtual momentum q of the gauge boson is spacelike, x is the familiar Bjorken variable. 



3. The Renormalisation Group and Asymptotic Freedom 

In this section we aim at providing a reasonably detailed introduction to the renormalisation group for- 
malism and the concept of running coupling which leads to the result that QCD has the property of 
asymptotic freedom. We start with a summary on how renormalisation works. 

In the simplest conceptual situation imagine that we implement regularisation of divergent inte- 
grals by introducing a dimensional cut-off A that respects gauge and Lorentz invariance. The dependence 
of renormalised quantities on A is eliminated by absorbing it into a redefinition of m (the quark mass: for 
simplicity we assume a single flavour here), the gauge coupling e (can be e in QED or e s in QCD) and the 

1/2 

wave function renormalisation factors Z q [ g for q and g, using suitable renormalisation conditions (that 
is precise definitions of m, g and Z). For example we can define the renormalised mass m as the position 



Fig. 6: 



of the pole in the quark propagator and, similarly, the normalisation Z q as the residue at the pole: 

Z q 

Propagator = — = — - — = + no — pole terms (20) 
pz — m z 

The renormalised coupling e can be defined in terms of a renormalised 3-point vertex at some specified 
values of the external momenta. We now become more specific by concentrating in the case of massless 
QCD. If we start from a vanishing mass at the classical (or "bare") level, mo = 0, the mass is not 
renormalised because it is protected by a symmetry, chiral symmetry. The conserved currents of chiral 
symmetry are axial currents: (77^75(7. The divergence of the axial current gives, by using the Dirac 
equation, 3^(^7^755) = Imq^q. So the axial current and corresponding axial charge are conserved 
in the massless limit. Since QCD is a vector theory we have not to worry about chiral anomalies in this 
respect. So one can choose a regularisation that preserves chiral symmetry besides gauge and Lorentz 
symmetry. Then the renormalised mass remains zero. The renormalised propagator has the form in 
eq.(^) with m = 0. 

The renormalised coupling e s can be defined from the renormalised 3-gluon vertex at a scale — fi 2 
(Fig. 7): 

V bare (p 2 ,q\r 2 ) = ZV ren (p 2 ,q 2 ,r 2 ), Z = Z Z J 2 Z V , e s = V ren (—fi 2 , — /i 2 , — /i 2 ) (21) 




r2 q 2 

Fig. 7: 

Here V\, are is what is obtained from computing the Feynman diagrams including, for example, the 
1-loop corrections at the lowest non trivial order (Vbare is defined so that it coincides with e s o in lowest 
order). It contains the cut-off A but does not know /x. Z is a factor that depends on the cut-off. Because 
of infrared singularities the defining scale fi cannot vanish. The negative value — /i 2 < is chosen to stay 



away from physical cuts (a gluon with negative virtual mass cannot decay). Similarly we can define Z g 
from the massless gluon propagator at the same scale — fi 2 (the vanishing mass of the gluon is guaranteed 
by gauge invariance). 

After computing all 1-loop diagrams in Fig. 7 we have: 

A 2 

Vbare(p 2 ,p 2 ,p 2 ) = ^Os [1 + Ca 0s ■ log + .. .] = 

A 2 -u 2 

= [l + ca,-log ? + ...]e s [l + ca s -log-£- + ...] = ZV ren (22) 

— fi z jr 

Note the replacement of eo with e in the second step. The definition of e s demands that one pre- 
cisely specifies what is included in Z. For this, in a given renormalisation scheme, a prescription is 
fixed to specify the finite terms that go into Z (i.e. the terms of order a s that accompany log A 2 ). 
Then V ren is specified and the renormalised coupling is defined from it according to eq.([2l|). For ex- 
ample, in the momentum subtraction scheme we define V ren (p 2 , p 2 , p 2 ) = e s + V& are (p 2 ,p 2 ,p 2 ) — 
Vbare{—^ 2 , — l^ 2 , — ^ 2 ), which is equivalent to say, at 1-loop, that all finite terms that do not vanish at 
p 2 = —fi 2 are included in Z. 

A crucial observation is that Vt, are depends on A but not on fi, which is only introduced when Z, 
V ren and hence a s are defined. (From here on, for shorthand, we write a to indicate either the QED 
coupling or the QCD coupling a s ). More in general for a generic Green function G, we similarly have: 



so that we have: 



or 



G bare (A 2 ,a ,p 2 ) = ZG ren (ii 2 ,a,p 2 ) (23) 

dGhare d 



d log fi 2 d log fi 2 



[ZG ren ] = (24) 



d da d 1 dZ 

1 d log /i 2 d log fi 2 da Z d log fi 
Finally the renormalisation group equation (RGE) can be written as: 



Z ^^-2 + 7^72^ + = (25) 



d log [i 2 da 

where 



+ (i{a)— + j(a)]Gren = (26) 



da 
d log /i 2 



/?(«) = iTTj (27) 



and 



7(«) = inr^2 (28) 



dlogZ 
d log ^i 2 

Note that (3(a) does not depend on which Green function G we are considering, but it is a property of 
the theory and the renormalisation scheme adopted, while 7(a) also depends on G. 

Assume that we want to apply the RGE to some hard process at a large scale Q, related to a Green 
function G that we can always take as adimensional (by multiplication by a suitable power of Q). Since 
the interesting dependence on Q will be logarithmic we introduce the variable t as : 

O 2 

t = log ^ (29) 
I 1 

Then we can write G ren = F(t, a, Xi) where x% are scaling variables (we often omit to write them in the 
following). In the naive scaling limit F should be independent of t. To find the actual dependence on t, 
we want to solve the RGE 

[—g- t + P( a )-fa + 7(a)] Gren = (30) 



with a given boundary condition at t = (or Q 2 = fi 2 ): F(0, a). 

We first solve the RGE in the simplest case that 7(0) = 0. This is not an unphysical case: for 
example, it applies to R e + e - where the vanishing of 7 is related to the non renormalisation of the electric 
charge in QCD (otherwise the proton and the electron charge would not exactly compensate: this will be 
better explained later). So we consider the equation: 

\-qI+ 0(a)—]G ren = (31) 

The solution is simply 

F(t,a) = F[0,a(t)] (32) 
where the "running coupling" a(t) is defined by: 

1 = / -Rf^ da (33) 
Ja p(pt!) 

Note that from this definition it follows that a(0) = a, so that the boundary condition is also satisfied. 
To prove that F[0, a(t)] is indeed the solution, we first take derivatives with respect of t and a (the two 
independent variables) of both sides of eq.(33). By taking d/dt we obtain 

1 = _J_^) (34 ) 
0(a(t) dt y ' 



We then take d/da and obtain 



1 1 da(t) 
= 1 — (35) 

/3(a) 0(a(t) da V ' 



These two relations make explicit the dependence of the running coupling on t and a 

da(t) 



dt 

da(t) _ (3{a{t)) 
da /3(a) 



(3(a(t)) (36) 

(37) 



Using these two equations one immediately checks that F[0, a(t)] is indeed the solution. 

Similarly, one finds that the solution of the more general equation with 7 / 0, eq.(30), is given 

by: 

fCt{t) ry( a ') 

F(t,a) = F[0,a(t)]exp / (38) 

Ja p{a') 

In fact the sum of the two derivatives acting on the factor F[0, a(t)] vanishes and the exponential is by 
itself a solution of the complete equation. Note that the boundary condition is also satisfied. 

The important point is the appearance of the running coupling that determines the asymptotic 
departures from scaling. The next step is to study the functional form of the running coupling. From 
eq. ([?(]) we see that the rate of change with t of the running coupling is determined by the (3 function. 
In turn /3(a) is determined by the \i dependence of the renormalised coupling through eq.(p7|). Clearly 
there is no dependence on (jl of the basic 3-gluon vertex in lowest order (order e). The dependence starts 
at 1-loop, that is at order e 3 (one extra gluon has to be emitted and reabsorbed). Thus we obtain that in 
perturbation theory: 

de n 



Recalling that a = e 2 /4tt, we have: 



a f a 2 oc 2e f e oc e 4 oc a 2 (40) 
a log ^ o log ;U Z 

Thus the behaviour of f3(a) in perturbation theory is as follows: 

(3(a) = ±ba 2 [l + b'a + ...] (41) 

Since the sign of the leading term is crucial in the following discussion, we stipulate that always b > 
and we make the sign explicit in front. By direct calculation at 1-loop one finds: 

QED : /3(a) ~ + ha 2 + b = V Nc(Q2)i (42) 

i 

where Nq = 3 for quarks and Nc = 1 for leptons and the sum runs over all fermions of charge Qe that 
are coupled. Also, one finds: 

QCD : [3(a) ba 2 + b = 11N c ~ 2llf (43) 
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where, as usual, nj is the number of coupled flavours of quarks (we assume here that nj < 16 so that 
b > in QCD). If a(t) is small we can compute (3(a(t)) in perturbation theory. The sign in front of b 
then decides the slope of the coupling: a(t) increases with t (or Q 2 ) if (3 is positive at small a (QED), or 
a(t) decreases with t (or Q 2 ) if (3 is negative at small a (QCD). A theory like QCD where the running 
coupling vanishes asymptotically at large Q 2 is called (ultraviolet) "asymptotically free". An important 
result that can be proven is that in 4 spacetime dimensions all and only non-abelian gauge theories can 
be asymptotically free. 

Going back to eq.(33) we replace /3(a) ~ ± ba 2 , do the integral and perform a simple algebra. 



We find 



and 



a 

QED : a t ~ - — (44) 

1 — bat 

QCD : a(t) ~ (45) 
1 + bat 

A slightly different form is often used in QCD. Defining 1/a = b log /x 2 / A.q CD we can write: 

111 

a(t) ~ --. = = = — (46) 

l + bt 6 i og ^ + Mog ^ blog^- 

QCD r QCD 

We see that a(t) decreases logarithmically with Q 2 and that one can introduce a dimensional parameter 
Aqcd that replaces /x. Often in the following we will simply write A for Aqcd, assuming that confusion 
with A = ultraviolet cut-off is avoided by the reader. Note that it is clear that A depends on the particular 
definition of a, not only on the defining scale fi but also on the renormalisation scheme (see, for example, 
the discussion in the next session). Through the parameter b, and in general through the j3 function, it 
also depends on the number nj of coupled flavours. It is very important to note that QED and QCD 
are theories with "decoupling": up to the scale Q only quarks with masses m << Q contribute to the 
running of a. This is clearly very important, given that all applications of perturbative QCD so far apply 
to energies below the top quark mass m t . For the validity of the decoupling theorem it is necessary that 
the theory where all the heavy particle internal lines are eliminated is still renormalisable and that the 
coupling constants do not vary with the mass. These requirements are true for the mass of heavy quarks 
in QED and QCD, but are not true in the electroweak theory where the elimination of the top would 
violate SU (2) symmetry (because the t and b left quarks are in a doublet) and the quark couplings to the 



Higgs multiplet (hence to the longitudinal gauge bosons) are proportional to the mass. In conclusion, 
in QED and QCD, quarks with m » Q do not contribute to rif in the coefficients of the relevant j3 
function. The effects of heavy quarks are power suppressed and can be taken separately into account. 
For example, in e + e~ annihilation for 2m c < Q < Im^ the relevant asymptotics is for nf = 4, while 
for 2m,5 < Q < 2m t nj = 5. Going accross the b threshold the function coefficients change, so the 
a(t) slope changes. But a(t) is continuous, so that A changes so as to keep constant a(t) at the matching 
point at Q ~ o(mb). The effect on A is large: approximately A5 ~ O.65A4. 



Note the presence of a pole in eqs.( |44[ , |45| ) at ±bat = 1, called the Landau pole, who realised 
its existence in QED already in the '50's. For [i ~ m e (in QED) the pole occurs beyond the Planck 
mass. In QCD the Landau pole is located for negative t or at Q < /x in the region of light hadron masses. 
Clearly the issue of the definition and the behaviour of the physical coupling in the region around the 
Landau pole is a problem that lies outside the domain of perturbative QCD. 

The non leading terms in the asymptotic behaviour of the running coupling can in principle be 
evaluated going back to eq.([4~ij) and computing b' at 2-loops and so on. But in general the perturba- 
tive coefficients of (3(a) depend on the definition of the renormalised coupling a (the renormalisation 
scheme), so one wonders whether it is worthwhile to do a complicated calculation to get b' if then it must 
be repeated for a different definition or scheme. In this respect it is interesting to remark that actually both 
b and b' are independent of the definition of a, while higher order coefficients do depend on that. Here 
is the simple proof. Two different perturbative definitions of a are related by a' ~ a(l + c\a + ...). 
Then we have: 



. r ua ua . 



da' da 
d log fi 2 d log fi 2 

±ba 2 (l + b'a + ...)(1 + 2 Cl a + ...) 
±6q' 2 (1 + b'a' + ... (47) 



which shows that, up to the first subleading order, j3{a') has the same form as /3(a). 
In QCD (Nq = 3) one has calculated: 
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2vr(33 - 2n f ) ( 

By taking b' into account one can write the expression of the running coupling at next to the leading 
order (NLO): 

a(Q 2 ) = a LO (Q 2 )[l - b'a LO (Q 2 )loglog^ + ...} (49) 

where a^Q = b log Q 2 /A 2 is the LO result. 

Summarizing, we started from massless classical QCD which is scale invariant. But we have seen 
that the procedure of quantisation, regularisation and renormalisation necessarily breaks scale invariance. 
In the quantum QCD theory there is a scale of energy, Aqcd, which from experiment is of the order of a 
few hundred MeV, its precise value depending on the definition, as we shall see in detail. Adimensional 
quantities depend on the energy scale through the running coupling which is a logarithmic function of 
Q 2 /A 2 . In QCD the running coupling decreases logarithmically at large Q 2 (asymptotic freedom), while 
in QED the coupling has the opposite behaviour. 



4. More on the Running Coupling 

In the previous section we have introduced the renormalised coupling a in terms of the 3-gluon vertex 
at p 2 = —fi 2 (momentum subtraction). The Ward identities of QCD then ensure that the coupling de- 
fined from other vertices like the qqg vertex are renormalised in the same way and the finite radiative 



corrections are related. But at present the universally adopted definition of a s is in terms of dimen- 
sional regularisation because of computational simplicity which is essential given the great complexity 
of present day calculations. So we now briefly review the principles of dimensional regularisation and the 
definition of Minimal Subtraction (MS) and Modified Minimal Subtraction (MS). The MS definition 
of a s is the one most commonly adopted in the literature and a value quoted for it is nomally referring 
to this definition. 

Dimensional Regularisation (DR) is a gauge and Lorentz invariant regularisation that consists 
in formulating the theory in D < 4 spacetime dimensions in order to make loop integrals ultraviolet 
finite. In DR one rewrites the theory in D dimensions (D is integer at the beginning, but then we will 
see that the expression of diagrams makes sense at all D except for isolated singularities). The metric 
tensor is extended into a D x D matrix = diag(l, — 1,— 1, 1) and 4- vectors are given by 
k^ = (k°, fc 1 , k ^ 1 ). The Dirac 7^ are f{D) x f(D) matrices and it is not important what is the 
precise form of the function f(D). It is sufficient to extend the usual algebra in a straightforward way 
like 7^7^ =-(£>- 2)Y or Tr( 7 ^) = f(D)g^. 

The physical dimensions of fields change in D dimensions and, as a consequence, the gauge cou- 
plings become dimensional eo = /^e, where e is adimensional, D = 4 — 2e and \i is a scale 
of mass (this is how a scale of mass is introduced in the DR of massless QCD!). The dimension of 
fields is determined by requiring that the action S = J dPxC is adimensional. By inserting for L 
terms like m^ty or m 2 (fy(j) or e^j^^A^ the dimensions of the fields and coupling are determined as: 
m, \P, (f), A^, e = 1, (D - 1) /2, (D - 2) /2, (D - 2) /2, (4 - D) /2, respectively. The formal expression 
of loop integrals can be written for any D. For example: 

dPk 1 _ r(2- D/2){-m 2 ) D l 2 - 2 



(2n) D (k 2 - m 2 ) 2 (4^/ 2 
For D = 4 — 2e one can expand using: 

r(e) = - - 1E + o(e), 7£ = 0.5772 (51) 

e 

For some Green function G, normalised to 1 in lowest order, (like V/e with V the 3-g vertex function at 
the symmetric point p 2 = q 2 = r 2 , considered in the previous section) we typically find at 1-loop: 

-u 2 1 

G bare = 1 + aoi-^f [c(- + lo g 4vr - lE ) + d + o(e)] (52) 
In MS one rewrites this at 1-loop accuracy (diagram by diagram: this is a virtue of the method): 

Z = 1 + a [c(- + log4vr - 7^)] 

_„ 2 

G ren = 1 + a [c log — 7T- + d] 

(53) 

In the original M S prescription only 1/e was subtracted (that clearly plays the role of a cutoff) and 
not also log47r and 7^. Later, since these constants always appear from the expansion of T functions 
it was decided to modify MS into MS. Note that the M S definition of a is different than that in the 
momentum subtraction scheme because the finite terms (those beyond logs) are different. In particular 
here 5G ren does not vanish at p 2 = —jjL 2 . 



The third coefficient of the QCD (3 function is also known in the M S prescription (recall that only 
the first two coefficients are scheme independent). Translated in numbers, for nj = 5 one obtains [Q|: 

/5(a) = -0.610a 2 [l + 1.261...- + 1.475.. .(-) 2 + ...] (54) 

7T 7T 



It is interesting to remark that the expansion coefficients are all of order 1, so that the MS expansion 
looks well behaved. 

5. Application to Hard processes 
5.1 R e + e - and Related Processes 

The simplest hard process is R e + e - that we have already started to discuss. R is adimensional and in 
perturbation theory is given by R = AfcEiQ^^sj.wheref = 1 + 0(05). We have already 
mentioned that for this process the "anomalous dimension" function vanishes: j(a s ) = because of 
electric charge non renormalisation by strong interactions. Let us review how this happens in detail. The 
diagrams that are relevant for charge renormalisation in QED at 1-loop are shown in Fig. 8. The Ward 
identity that follows from gauge invariance in QED imposes that the vertex (Zy) and the self-energy 
(Zf) renormalisation factors cancel and the only divergence remains in Z 7 , the vacuum polarization of 
the photon. So the charge is only renormalised by the photon blob, hence it is universal (the same factor 
for all fermions, independent of their charge) and is not affected by QCD at 1-loop. It is true that at 
higher orders the photon vacuum polarization diagram is affected by QCD (for example, at 2-loops we 
can exchange a gluon between the quarks in the photon loop) but the renormalisation induced by the 
vacuum polarisation diagram remains independent of the nature of the fermion to which the photon line 
is attached. The gluon contributions to the vertex {Zy) and to the self-energy (Zf) cancel because they 
have exactly the same structure as in QED, and there is no gluon contribution to the lowest order photon 
vacuum polarisation blob. So j(a s ) = 0. 





At 1-loop the diagrams relevant for the computation of R are shown in Fig. 9. There are virtual 
diagrams and real diagrams with one additional gluon in the final state. Infrared divergences cancel 
between the interference term of the virtual diagrams and the absolute square of the real diagrams, 
according to the Bloch-Nordsieck theorem. Similarly there are no mass singularities, in agreement with 
the Kinoshita-Lee-Nauenberg theorem, because the initial state is purely leptonic and all degenerate 
states that can appear at the given order are included in the final state. Given that j(a s ) = the RGE 
prediction is simply given, as we have already seen, by F(t, a s ) = F[0, a s (t)]. This means that if we 
do, for example, a 2-loop calculation, we must obtain a result of the form: 

F(t,a s ) = 1 + cia s (l — ba s t) + C20? s + o(a^) (55) 

In fact we see that this form, taking into account that from eq.(p3|) we have: 

Qt 

a s (t) ~ i + b a t ~ ~ bast + ""^ ^ 



Fig. 9: 



can be rewritten as 

F(t,a s ) = 1 + cia a (t) + c 2 a 2 s {t) + o(a 3 s (t)) = F[0,a s {t)} (57) 

The content of the RGE prediction is, at this order, that there are no a s t and (a s t) 2 terms (the leading 
log sequence must be absent) and the term of order a 2 t has the coefficient that allows to reabsorb it in 
the transformation of a s into a s (t). 

At present the first 3 coefficients have been computed in the MS scheme [[J. Clearly c\ = 1/tt 
does not depend on the definition of a s but C2 and C3 do. The subleading coefficients also depend on 
the scale choice: if instead of expanding in a s (Q) we decide to choose a s (Q/2) the coefficients C2 
and C3 change. In the MS scheme, for 7-exchange and nj = 5, which are good approximations for 

2mf, « Q « mz, one has: 

F[0,a s (t)} = 1 + ^ + 1.409... (^) 2 - 12.8...(^) 3 + ... (58) 

TT TT TT 

Similar perturbative results at 3-loop accuracy also exist for Rz = T(Z — ► hadrons) /T(Z — > 
leptons), R T = T(t — > v T + hadrons) /T(t — > v T + leptons), etc. We will discuss these results later 
when we deal with measurements of a s . 

The perturbative expansion in powers of a s (t) takes into account all contributions that are sup- 
pressed by powers of logarithms of the large scale Q 2 ("leading twist" terms). In addition there are 
corrections suppressed by powers of the large scale Q 2 ("higher twist" terms). The pattern of power cor- 
rections is controlled by the light-cone Operator Product Expansion (OPE) which (schematically) leads 
to: 

77 f ^ m 2 <Q|Tr[F^F^]|0> < 0|O 6 |0 > , 

F = pert. + r 2 — 7j + r 4 — L -^ I J h ••• + r 6 >— J h ... (59) 



Here m 2 generically indicates mass corrections, notably from b quarks, for example (t quark mass cor- 
rections only arise from loops, vanish in the limit mt — ► 00 and are included in the coefficients as those in 
eq.(^) and the analogous ones for higher twist terms), = J2a F^i,t A , 0§ is typically a 4-fermion 
operator, etc. For each possible gauge invariant operator the corresponding power of Q 2 is fixed by 
dimensions. 

We now consider the light-cone OPE in some more detail. R e + e - ~ n(Q 2 ) where H(Q 2 ) is the 
scalar spectral function related to the hadronic contribution to the imaginary part of the photon vacuum 
polarization T^ u : 

T^ u = (-g^uQ 2 + qnQv)n(Q 2 ) = J expiqx < 0\Jl(x)J u (0)\0 > dx = 



= £< 0|jj(0)|n X n| J„(0)|0 > (27r) 4 5 4 (g - p n ) 

n 



(60) 



For Q 2 —> oo the x 2 — > region is dominant. To all orders in perturbation theory the OPE can be proven. 
Schematically, dropping Lorentz indices, for simplicity, near x 2 ~ we have: 

oo 

4(x)J u (0) = I(x 2 ) + E(x 2 )Y / c n (x 2 )x^...x^- O^JO) + 

n=0 

+ less sing, terms (61) 



Here I(x 2 ), E(x 2 ),..., c n (x 2 ) are c-number singular functions, O n is a string of local operators. E(x ) 
is the singularity of free field theory, I(x 2 ) and c n (x 2 ) contain powers of log \ix in interaction. Some O n 
are already present in free field theory, other ones appear when interactions are switched on. n(Q 2 ) is 
related to the Fourier transform. Less singular terms in x 2 lead to power suppressed terms in l/Q 2 . The 
perturbative terms come from I{x 2 ) which is the leading twist term. The logarithmic scaling violations 
induced by the running coupling are the logs in I(x 2 ). 

5.2 The Final State in e+e Annihilation 

Experiments on e + e~ annihilation at high energy provide a remarkable possibility of systematically test- 
ing the distinct signatures predicted by QCD for the structure of the final state averaged over a large 
number of events. Typical of asymptotic freedom is the hierarchy of configurations emerging as a conse- 
quence of the smallness of a s {Q 2 ). When all corrections of order a s {Q 2 ) are neglected one recovers the 
naive parton model prediction for the final state: almost collinear events with two back-to-back jets with 
limited transverse momentum and an angular distribution as (1 + cos 2 6) with respect to the beam axis 
(typical of spin 1/2 parton quarks: scalar quarks would lead to a sin 2 9 distribution). At order a s (Q 2 ) 
a tail of events is predicted to appear with large transverse momentum p? ~ Q/2 with respect to the 
thrust axis (the axis that maximizes the sum of the absolute values of the longitudinal momenta of the 
final state particles). This small fraction of events with large px mostly consists of three-jet events with 
an almost planar topology. The skeleton of a three-jet event, at leading order in a s (Q 2 ), is formed by 
three hard partons qqg, the third being a gluon emitted by a quark or antiquark line. The distribution of 
three-jet events is given by: 

1 da = 2a s x\ + x\ 

a dx\dx2 3ir (1 — X\)(l — x^) 

here x\^ refer to energy fractions of massless quarks: Xi = 2E; L j ^fs with x\ + X2 + X3 = 2. At order 
a 2 (Q 2 ) a hard perturbative non planar component starts to build up and a small fraction of four-jet events 
qqgg or qqqq appear, and so on. 

For precise testing and for measuring a s a quantitatively specified definition of jet counting must 
be introduced which must be infrared safe (i.e. not altered by soft particle emission or collinear splittings 
of massless particles) in order to be computable at parton level and as much as possible insensitive to the 
transformation of partons into hadrons. One introduces a resolution parameter y cut and a suitable pair 
variable, for example: 

min(E?,E?)(l -cosdjj) 

Vij = (63) 

s 

The particles i,j belong to different jets for > y cut . Clearly the number of jets becomes a function 
of y C ut- there are more jets for smaller y cu t. Measurements of a s (Q 2 ) have been performed starting 
from jet multiplicities, the largest error coming from the necessity of correcting for non-perturbative 
hadronisation effects. 



5.3 Deep Inelastic Scattering 

Deep Inelastic Scattering (DIS) processes have played and still play a very important role for our under- 
standing of QCD and of nucleon structure. This set of processes actually provides us with a rich labora- 
tory for theory and experiment. There are several structure functions that can be studied, Fi(x, Q 2 ), each 
a function of two variables. This is true separately for different beams and targets and different polariza- 
tions. Depending on the charges of 1 and 1' we can have neutral currents (7,Z) or charged currents in the 
l'-l channel (Fig. 6). In the past DIS processes were crucial for establishing quarks and gluons as partons 
and QCD as the theory of strong interactions. At present DIS is very important for quantitative studies 
and tests of QCD. The theory of scaling violations for totally inclusive DIS structure functions, based on 
operator expansions and renormalization group techniques, is crystal clear and the predicted Q 2 depen- 
dence can be tested at each value of x. The measurement of quark and gluon densities in the nucleon, 
as functions of x at some reference value of Q 2 , which is an essential starting point for the calculation 
of all relevant hadronic hard processes, is performed in DIS processes. At the same time one measures 
a s (Q 2 ) and the DIS values can be compared with those obtained from other processes. At all times new 
theoretical challenges arise from the study of DIS processes. Recent examples are the so-called "spin 
crisis" in polarized DIS and the behaviour of singlet structure functions at small x as revealed by HERA 
data. In the following we will review the past successes and the present open problems in the physics of 
DIS. 

The cross-section a ~ L^ V W^ U is given in terms of the product of a leptonic (L^ v ) and a hadronic 
(W^) tensor. While W is simple and easily obtained from the lowest order electroweak (e-w) vertex 
plus QED radiative corrections, the complicated strong interaction dynamics is contained in W^ u . The 
latter is proportional to the Fourier transform of the forward matrix element between the nucleon target 
states of the product of two e-w currents: 

= j dx expiqx < p\ J„(0)|p > (64) 

Structure functions are defined starting from the general form of W^ v given Lorentz invariance and 
current conservation. For example, for e-w currents between unpolarized nucleons we have: 

w i I^Mxr / , i mv mv W 2 {v,Q 2 ) 



q* q* ' ' q* m 2 



^^uX P p A q p W^Q 2 ) 



W 3 is absent for pure vector currents. In the limit Q 2 » m 2 , x fixed, the structure functions obey 
approximate Bjorken scaling which in reality is broken by logarithmic corrections that can be computed 
in QCD: 

mW x (y,Q 2 ) -► Fi(ar) 
vW 2 , 3 (v,Q 2 )^F 2>3 (x) (65) 

The 7 — N cross-section is given by (Wj = Wi(Q 2 ,u)): 

^ = ^•[23^/2^ + 00^/2^] (66) 
while for the v — N or v — N cross-section one has: 

iSr = ^ U^ 2 ) 2 -[2sm 2 e/2 Wl + cos 2 9/2W 2 ±^±^ sin 2 9/2W 3 ] (67) 
(Wi for photons, v and v are all different, as we shall see in a moment). 



In the scaling limit the longitudinal and transverse cross sections are given by: 




0"T = ORH + &LH 



(68) 



where L, RH, LH refer to the helicity 0, 1, -1, respectively, of the exchanged gauge vector boson. 

In the '60's the demise of hadrons from the status of fundamental particles to that of bound states 
of constituent quarks was the breakthrough that made possible the construction of a renormalisable field 
theory for strong interactions. The presence of an unlimited number of hadrons species, many of them 
with large spin values, presented an obvious dead-end for a manageable field theory. The evidence for 
constituent quarks emerged clearly from the systematics of hadron spectroscopy. The complications 
of the hadron spectrum could be explained in terms of the quantum numbers of spin 1/2, fractionally 
charged, u, d and s quarks. The notion of colour was introduced to reconcile the observed spectrum 
with Fermi statistics. But confinement that forbids the observation of free quarks was a clear obstacle 
towards the acceptance of quarks as real constituents and not just as fictitious entities describing some 
mathematical pattern (a doubt expressed even by Gell-Mann at the time). The early measurements at 
SLAC of DIS dissipated all doubts: the observation of Bjorken scaling and the success of the naive (not 
so much after all) parton model of Feynman imposed quarks as the basic fields for describing the nucleon 
structure (parton quarks). 

In the language of Bjorken and Feynman the virtual 7 (or, in general, any gauge boson) sees the 
quark partons inside the nucleon target as quasi-free, because the (Lorentz dilated) QCD interaction 
time is much longer than r 7 ~ l/Q- Since the virtual photon 4-momentum is spacelike, we can go to 
a Lorentz frame where E 1 = (Breit frame). In this frame q = (E 1 = 0; 0, 0, Q) and the nucleon 
momentum, neglecting the mass m << Q, is p = (Q/2x; 0, 0, —Q/2x) (note that this correctly gives 
x = Q 2 /2(p ■ q)). Consider (Fig. 10) the interaction of the photon with a quark carrying a fraction y 
of the nucleon 4-momentum: p g = yp (we are neglecting the transverse components of p g which are of 
order m). The incoming parton with p q = yp absorbs the photon and the final parton has 4-momentum 
p' q . Since in the Breit frame the photon carries no energy but only a longitudinal momentum Q, the 
photon can only be absorbed by those partons with y = x: then the longitudinal component of p q = yp 
is —yQ/2x = —Q/2 and can be flipped into +Q/2 by the photon. As a result, the photon longitudinal 
momentum +Q disappears, the parton quark momentum changes of sign from —Q/2 into +Q/2 and the 
energy is not changed. So the structure functions are proportional to the density of partons with fraction x 
of the nucleon momentum, weighted with the squared charge. Also, recall that the helicity of a massless 
quark is conserved in a vector (or axial vector) interaction. So when the momentum is reversed also 
the spin must flip. Since the process is collinear there is no orbital contribution and only a photon with 
helicity ±1 (transverse photon) can be absorbed. If partons were spin zero only longitudinal photons 
would instead contribute. 



-Q/2 




Fig. 10: 



Using these results, which are maintained in QCD at leading order, the quantum numbers of the 
quarks were confirmed by early experiments. The observation that R = ai/cr —> implies that the 



charged partons have spin 1/2. The quark charges were derived from the data on the electron and neutrino 
structure functions: 



F ep = A/9u(x) + l/9d(x) + ; F en = 4/9d(x) + l/9u{x) + .... 

Fup = F Pn = 2d{x) + ; F vn = F Vp = 2u(x) + (69) 

where F ~ 2F\ ~ F2/X and u(x), d{x) are the parton number densities in the proton (with fraction x of 
the proton longitudinal momentum), which, in the scaling limit, do not depend on Q 2 . The normalisation 
of the structure functions and the parton densities are such that the charge relations hold: 

f\u(x) -u(x)]dx = 2, f [d{x) - d(x)]dx = 1, / [s{x] - s{x))dx = (70) 
Jo Jo Jo 

Also it was proven by experiment that at values of Q 2 of a few GeV 2 , in the scaling region, about half 
of the nucleon momentum, given by the momentum sum rule: 

I E(gt(aO+ft(aO) + 9{x)]xdx = 1 (71) 
Jo ; 

is carried by neutral partons (gluons). 

In QCD there are calculable log scaling violations induced by a s (t). The parton rules just intro- 
duced can be summarised in the formula: 

F(x,t) = f'dy^-a point (-,a s (t)) + o(-L) (72) 
Jx y y Q z 

Before QCD corrections (T po int = e 2 5(x/y — 1) and F = e 2 qo(x) (here we denote by e the charge of 
the quark in units of the positron charge, i.e. e = 2/3 for the u quark). QCD modifies a po i n t at order 
a s via the diagrams of Fig. 1 1 . Note that the integral is from x to 1 , because the energy can only be lost 
by radiation before interacting with the photon (which eventually wants to find a fraction x, as we have 
explained). From a direct computation of the diagrams one obtains a result of the following form: 

a point {z,a s (t)) ~ e 2 [5(z-l) + ^(t-P(z) + /(*))] (73) 



For y > x the correction arises from diagrams with real gluon emission. Only the sum of the 
two real diagrams in Fig. 1 1 is gauge invariant, so that the contribution of one given diagram is gauge 
dependent. There is a special form of axial gauge, called physical gauge, where, among real diagrams, 
the diagram of Fig. 11(c) gives the whole t-proportional term. It is obviously not essential to go to this 
gauge, but this diagram has a direct physical interpretation: a quark in the proton has a fraction y > x of 
the parent 4-momentum; it then radiates a gluon and looses energy down to a fraction x before interacting 
with the photon. The log arises from the virtual quark propagator, according to the discussion of collinear 
mass singularities in eq.(17). In fact in the massless limit one has: 

11 -11 



propagator 



r 2 (k-h) 2 2E k E h l-cosfl 
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" AE^ ' sin 2 9/2 ^ & ( } 

where px is the transverse momentum of the virtual quark. So the square of the propagator goes like 
1/Pt- But there is ap^ factor in the numerator, because in the collinear limit, when 9 = and the initial 
and final quarks and the emitted gluon are all aligned, the quark helicity cannot flip (vector interaction) 
and the real gluon cannot have but ±1 helicity. So the cross-section behaves as: 

a ~ / — dp 2 T ~ logQ 2 (75) 
J p T 




Actually the log should be read as log Q 2 /m 2 because in the massless limit a genuine mass singularity 
appears. In fact the mass singularity connected with the initial quark line is not cancelled because we 
do not have the sum of all degenerate initial states, but only a single quark. But in correspondence to 
the initial quark we have the (bare) quark density qo(y) that appear in the convolution integral. This is a 
non perturbative quantity that is determined by the nucleon wave function. So we can factorize the mass 
singularity in a redefinition of the quark density: we replace qo(y) — > q(y, t) = qo(y) + Ag(y, t) with: 



A*,.) = ?«f#^.pfi (76) 
2vr J x y y 

Here the factor of t is a bit symbolic: it stands for log Q 2 /km 2 and what we exactly put below Q 2 
depends on the definition of the renormalised quark density, that also determines the exact form of the 



finite term f(z) in eq.(73). 



The effective parton density q(y,t) that we have defined is now scale dependent. In terms of this 
scale dependent density we have the following relations, where we have also replaced the fixed coupling 
with the running coupling according to the prescription derived from the RGE: 

F(x,t) = f 1 dy q J^le 2 [5(- - 1) + ^/(-))] = e 2 q(x,t) + o(a s (t)) 
Jx y y 2vr y 

d a s (t) f 1 q(y,t) ,x 2 

-rrtfOM) = -= — / dy P(-) + o{a s {t) ) (77) 

at 2n J x y y 



We see that in lowest order we reproduce the naive parton model formulae for the structure functions in 
terms of effective parton densities that are scale dependent. The evolution equations for the parton densi- 
ties are written down in terms of kernels (the "splitting functions") that can be expanded in powers of the 
running coupling. At leading order, we can interpret the evolution equation by saying that the variation 
of the quark density at x is given by the convolution of the quark density at y times the probability of 
emitting a gluon with fraction x/y of the quark momentum. 

It is interesting that the integro-differential QCD evolution equation for densities can be trans- 
formed into an ordinary differential equation for Mellin moments. The moment /„ of a density f(x) is 
defined as: 

fn = ^ dxx^fix) (78) 
JO 



By taking moments of both sides of the second of eqs.(77) one finds, with a simple interchange of the 
integration order, the simpler equation for the nth moment: 



j t qn{t) = -f^ ■ P n ■ ?„(*) (79) 



To solve this equation we observe that: 



, qn(t) Pn /"* , s , Pn f as(t) da' 



g„(0) 27r Jo 27r J as -ba 

where we used eq.(|3^) to change the integration variable from dt to da(t) (denoted as da') and 13(a) ~ 
—ba 2 + .... Finally the solution is: 

?»(*) = ["^]^-«n(0) (81) 

The connection of these results with the RGE general formalism occurs via the light cone OPE 
(recall eq.(64) for and eq.([H]) for the OPE of two currents). In the case of DIS the c-number term 

I(x 2 ) does not contribute, because we are interested in the connected part < p|...|j> > — < 1 1 >. 

The relevant terms are: 

oo 

jJ(s)J„(0) = E{x 2 ) c n (x 2 )x^ 1 ...x lln ■ 0^.^(0) + less sing, terms (82) 

n=0 

A formally intricate but conceptually simple argument (Ref.[^jJ|, page 28) based on the analiticity prop- 
erties of the forward virtual Compton amplitude shows that the Mellin moments M n of structure func- 
tions are related to the individual terms in the OPE, precisely to the Fourier transform c n (Q 2 ) (we will 
write it as c n (t, a)) of the coefficient c n {x 2 ) times a reduced matrix element h n from the operators O n : 

< p|ojJi.../i n (o)b >= Kp^-p^- 



c n < p\O n \p >^ M n = f 1 dxx n ~ l F(x) 

Jo 



(83) 



Since the matrix element of the products of currents satisfy the RGE so do the moments M n . Hence the 
general form of the Q 2 dependence is given by the RGE solution (see eq.(|3^)): 



M n (t,a) = c n [0,a(t)]exp / J^±da' ■ h n (a) (84) 



In lowest order, identifying in the simplest case M n with q n , we have: 



Inipc) = ^a + (3(a) = -ba 2 + ... (85) 

Z7T 



and 



Qn(t) = q n (0) exp f ^ = ■ q n (0) (86) 

Ja pa a s (t) 



<■«(*) ln (a') , _as 
(3(a>) a l a s (ty 
which exactly coincides with eq.(81). 

Up to this point we have implicitely restricted our attention to non-singlet (under the flavour group) 
structure functions. The Q 2 evolution equations become non diagonal as soon as we take into account 
the presence of gluons in the target. In fact the quark which is seen by the photon can be generated by a 
gluon in the target (Fig. 12). 



The quark evolution equation becomes: 

*(x t) - ^ 

dt qA ' > ~ 2vr "'O-™ ' 2vr 




Fig. 12: 



where we introduced the shorthand notation: 



1 dy^ . p(E) 
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(88) 



(it is easy to check that the convolution, like an ordinary product, is commutative). At leading order, 
the interpretation of eq.(|87|) is simply that the variation of the quark density is due to the convolution of 
the quark density at a higher energy times the probability of finding a quark in a quark (with the right 
energy fraction) plus the gluon density at a higher energy times the probability of finding a quark (of 
the given flavour i) in a gluon. The evolution equation for the gluon density, needed to close the system, 
can be obtained by suitably extending the same line of reasoning to a gedanken probe sensitive to colour 
charges, for example a virtual gluon. The resulting equation is of the form: 
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(89) 



The explicit form of the splitting functions in lowest order can be explicitly derived from the QCD 
vertices. They are a property of the theory and do not depend on the particular process the parton density 
is taking part into. The results are: 

„2 



P 



HI 



4 r 1 + x z 3 . 

-h ; 1- -<*(1 - z) + o(a s 

3 L (l-x+ 2 v n v 



P, 



41 + (l-x) 2 
3 x 



o(a s ) 



P 



1 



Q9 



[x 2 + (1 - x) 2 ] + o(a s ) 



99 



6f 



+ 



+ x(l-z)] + 



33-2nf. , . , 
L S(1 - x) + o(a s ) 



(1 — x) + x ' s 6 

The "+" distribution is defined as, for a generic non singular weight function f(x): 



fix) 



o (1 - x)+ 
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dx 
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(91) 



The 5(1 — x) terms arise from the virtual corrections to the tree diagram. Their coefficient can be simply 
obtained by imposing the validity of charge and momentum sum rules. In fact, from the request that the 
charge sum rules in eq.(f70|) are not affected by the Q 2 dependence one derives that 

■i 



P qq (x)dx 







(92) 



which can be used to fix the coefficient of the 5(1 — x) terms of P qq . Similarly, by taking the t-derivative 
of the momentum sum rule in eq.([7j]) and imposing its vanishing for generic q-i and g, one obtains: 

f\p qq (x) + P gq (x)]xdx = 0, [ 1 {2n f P qg {x) + P gg (x)]xdx = 0, (93) 
Jo Jo 

At higher orders the evolution equations are easily generalised but the calculation of the splitting 
functions rapidly becomes very complicated. The splitting functions are completely known at NLO 
accuracy: a s P ~ a s P\ + a 2 s P2 + .... More recently the NNLO results P3 have been derived in 
analytic form for the first few moments [f7|]. The full NNLO calculation is in progress and could be 
finished soon. 

The scaling violations are clearly observed by experiment and their pattern is very well reproduced 
by QCD fits at NLO. Examples are seen in Figs. 13(a-d) These fits provide an impressive confir- 
mation of a quantitative QCD prediction, a measurement of qi(x, Qq) and g(x, Qq) at some reference 
value Qq of Q 2 and a precise measurement of a s (m 2 z ). At small x and large but fixed Q 2 when terms of 
order (a s log l/x) n become of order 1 and cannot be neglected, the validity of the NLO or NNLO ap- 
proximations should breakdown. However, the small x data collected by HERA can be fitted reasonably 
even at the smallest measured values of x by the NLO QCD evolution equations, so that there is no neat 
evidence in the data for departures. But the extracted gluon density (which is dominant at small x) and 
the fitted value of a s (m? z ) could be biassed by this effect if data at too small values of x are included. 



5.4 Factorisation and the QCD Improved Parton Model 

The parton densities defined and measured in DIS are used to compute hard processes initiated by 
hadronic collisions via the Factorisation Theorem (FT). Suppose you have a hadronic process of the 
form hi + hi — ► X + all where X is some triggering particle or pair of particles which specify the 
large scale Q 2 relevant for the process, in general somewhat but not much smaller than s, the total c.o.m. 
squared mass. For example, in pp or pp collisions, X can be a W or a Z or a virtual photon with 
large Q 2 , or a jet at large transverse momentum pr, or a pair of heavy quark-antiquark of mass M. By 
"all" we mean a totally inclusive collection of gluons and light quark pairs. The FT implies that for the 
cross-section or some other sufficiently inclusive distribution we can write the expression: 

<t(s,t) = / dx 1 dx 2 piA(xi,Q 2 )p2B(x2,Q 2 )<7AB(xix 2 s,T) (94) 
AB J 

Here r = Q 2 /s is a scaling variable, pic are generic parton-C densities inside the hadron hi, gab is the 
partonic cross-section for parton-A + parton-B^ X + all'. This result is based on the fact that the mass 
singularities that are associated with the initial legs are of universal nature, so that one can reproduce 
the same modified parton densities, by absorbing these singularities into the bare parton densities, as 
in deep inelastic scattering. Once the parton densities and a s are known from other measurements the 
prediction of the rate for a given hard process is obtained with no free parameters. The NLO calculation 
of the reduced partonic cross-section is needed in order to correctly specify the scale and in general the 
definition of the parton densities and of the running coupling in the leading term. The residual scale and 
scheme dependence is often the most important source of theoretical error. In the following we consider 
a few examples. 

A comparison of data and predictions on the production of jets at large and px in pp or pp 
collisions is shown in Fig. 14 [Q]. 

This is a particularly significant test because the rates at different c.o.m. energies and, for each 
energy, at different values of px span over many orders of magnitude. Also the corresponding values of 
i/i and px are large enough to be well inside the perturbative region. The overall agreement of the data 
from ISR, UA1,2 and CDF and DO is spectacular. Only at very large px there might be some problem 
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Fig. 13: A recent NLO fit of scaling violations from Ref. [ph, for different x ranges, as functions of Q 2 




according to CDF and, to a lesser extent, to DO. A harder gluon at large x and inclusion of systematic 
errors can alleviate or eliminate the problem. This issue will be clarified in the near future by the Runll 
results at the Tevatron. Similar results also hold for the production of photons at large pr. The collider 
data, shown in Fig. 15 [Q], are in fair agreement with the theoretical predictions. A less clear situation is 
found with fixed target data. Here, first of all, the experimental results show some internal discrepancies. 
Also, the pt accessible values being smaller, the theoretical uncertainties are larger. But it is true that 
the agreement is poor, so that the necessity of an "intrinsic" transverse momentum of partons inside the 
hadron of over 1 GeV has been claimed, which I do not think one is free to introduce. 

For heavy quark production at colliders the agreement is very good for top production at the 
Tevatron (Fig. 16). To present a rare example of a puzzle, the bottom production at the Tevatron remains 
problematic. The total rate and the pt distribution of b quarks observed at CDF is in excess of the 
prediction, up to the largest measured values of pr (fig. 17). It is true that this is a complicated case, with 
different scales being present at the same time: y^i, pr, m^. Also some non-perturbative ingredients, 
like fragmentation functions, are present in the calculation and could be in part reponsible for the effect. 
But no clear explanation of this phenomenon is available. Probably it occurs through an unfavourable 
conspiracy of several small effects. 

Drell-Yan processes, including lepton pair production via virtual 7, W or Z exchange, offer a good 
opportunity to test QCD. The process is quadratic in parton densities, and the final state is totally inclu- 
sive, while the large scale is specified and measured by the invariant mass squared Q 2 of the lepton pair 
which itself is not strongly interacting. The QCD improved parton model leads directly to a prediction 
for the total rate as a function of Q 2 . The value of the LO cross-section is inversely proportional to the 
number of colours Nc because a quark of given colour can only annihilate with an antiquark of the same 
colour to produce a colourless lepton pair. The order a s (Q 2 ) corrections to the total rate were computed 
long ago and found to be particularly large, when the quark densities are defined from the structure func- 
tion F2 measured in DIS at q 2 = —Q 2 . The ratio a corr /aLO of the corrected and the Born cross-sections, 
was called K-factor (by me), because it is almost a constant in rapidity. More recently also the NLO full 
calculation of the K-factor was completed, a very remarkable calculation [§]. The QCD predictions can 
be best tested for W and Z production at CERN and Tevatron energies. Q ~ rnw,z is large enough to 
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Fig. 15: Single photon production in pp colliders as function of pr 



make the prediction reliable (a not too large K-factor) and the ratio y/r = Q/y/s is not too small. Recall 
that in lowest order x\X2S = Q 2 so that the parton densities are probed at x values around y/r. We have 
y/r = 0.13 - 0.15 at y/1 = 630 GeV (CERN SppS Collider) and y/r = 0.041 - 0.052 at the Tevatron. 
In this respect the prediction is more delicate at the LHC, where y/r ~ 5.9 — 6.5 • 10~ 3 . One comparison 
of the experimental total rates at the Tevatron with the QCD predictions is shown in Fig. 18, together 
with the expected rates at the LHC (based on the structure functions obtained in H). 

The calculation of the W/Z px distribution has been a classic problem in QCD. For large pr, for 
example px ~ o{mw), the px distribution can be reliably computed in perturbation theory, which was 
done up to NLO in the late '70's and early '80's. A problem arises in the intermediate range Aqcd << 
Pt « mw, where the bulk of the data is concentrated, because terms of order a s (p T )logm 2 v /p T 
become of order 1 and should included to all orders. At order a s we have: 

' da{) (1 + A)5(pl) + ^log™f +r^+ D(pl) (95) 



OQdp 2 T p\ p\ + (p^ J+ 

where A, B, C, D are coefficients of order a s . The "+" distribution is defined in complete analogy with 
eq.©: 

2 2 
PTMAX , , , , fPTMAX , 



g(z)f(z) + dz = / [g(z)-g(0)]f(z)dz (96) 

o Jo 

The content of this at first sight mysterious definition is that the singular "+" terms do not contribute to 

the total cross-section. In fact for the cross-section the weight function g{z) = 1 and we obtain: 

o = a [(l + A) + / D(z)dz] (97) 

jo 

The singular terms, of infrared origin, are present at the non completely inclusive level but disappear in 
the total cross-section. These singularities are proven to exponentiate and lead to the following expres- 
sion: 

1 duo f d 2 b 
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f — exp(-ib-p T )(l + A)expS(b) (98) 
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Fig. 17: The b production pr distribution at the Tevatron pp collider (CDF collaboration). The cross-section for pr > p™™ is 
plotted as function of p™" 1 . 



with: 

S{b)= ^ [e xp^-6-l][^log^ + ^] (99) 

At large pj- the LO perturbative expansion is recovered. At intermediate p? the infrared pj- singularities 
are resummed (the Sudakov log terms, which are typical of vector gluons, are related to the fact that for 
a charged particle in acceleration it is impossible not to radiate, so that the amplitude for no soft gluon 
emission is exponentially suppressed). However this formula has problems at small pr, for example, 
because of the presence of a s under the integral for S(b): presumably the relevant scale is of order k\. 
So it must be completed by some non perturbative ansatz or an extrapolation into the soft region. All 
the exercise has been extended to NLO accuracy, where one starts from the perturbative expansion at 
order a 2 s , and generalises the resummation to also include NLO terms of order a s (pj,) 2 logm^/pj,. 
The comparison with the data is very impressive. In Fig. 19 we see the px distribution as predicted in 
QCD (with a number of variants that mainly differ in the approach to the soft region) compared with 
some recent data from the DO Collaboration at the Tevatron [jlOl]. 
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Fig. 18: Data vs. theory for W and Z production at the Tevatron (y/s = 1.8 TeV) together with the corresponding predictions 
for the LHC (y/s= 1.4 TeV) 



6. Measurements of a s 

The most precise and reliable measurements of a s (m z ) in the MS definition are obtained from e + e~ 
colliders (in particular LEP) and from Deep Inelastic Scattering. 

6.1 a s from e + e colliders 

The main methods at e + e~ colliders are: a) Inclusive hadronic Z decay, Ri, ah, Yz- b) Inclusive hadronic 
r decay, c) Event shapes and jet rates. 

As we have seen, for a quantity like Ri we can write a general expression of the form: 

Ri = z ;^ hadrons l „ r^(i + s QCD + 5 NP ) + ... (ioo) 

Z, t — > leptons) 

where R EW is the electroweak-corrected Born approximation, 5qcd, <5/vp are the perturbative (logarith- 
mic) and non perturbative (power suppressed) QCD corrections. If we consider measurement at the Z, 
from Ri only, assuming the standard electroweak theory, one finds [jll|] : 

a s {m z ) = 0.123 ±0.004 (101) 

Better, one can use all info from Ri, Yz = 3T; + Yh and ah = 12irYiYh/ {rn z Y 2 z ). From these observ- 
ables together one obtains: 

a s (m z ) = 0.120 ±0.003 (102) 
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Fig. 19: QCD predictions for the W pr distribution compared with recent DO data at the Tevatron (y/s = 1.8 TeV) 



By adding all other electroweak precision electroweak tests (in particular mw) one finds: 

a s (m z ) = 0.118 ±0.003 (103) 

The final error is predominantly theoretical and is dominated by our ignorance on m# and from higher 
orders in the QCD expansion. 

We now consider the measurement of a s (mz) from r decay. R T has a number of advantages 
that, at least in part, tend to compensate for the smallness of m T = 1.777 GeV. First, R T is maximally 
inclusive, more than R e + e -(s), because one also integrates over all values of the invariant hadronic 
squared mass: 

1 f m r d,S , S N n 

R r = - -^(1 - -^) 2 ImU T (s) (104) 
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Analyticity can be used to transform the integral into one on the circle at \s\ = m 2 

2m J\ 3 \= m % m% m 2 T 

Also, the factor (1 — ^r) 2 is important to kill the sensitivity the region Res = m 2 . where the physical cut 
and the associated thresholds are located. Still the quoted result (a combination of ALEPH and OPAL 
analyses) looks a bit too precise: 

a s (m z ) = 0.1181 ± 0.0007(exp) ± 0.003(t/i) (106) 

This precision is obtained by taking for granted that corrections suppressed by 1/ m 2 are negligible. This 
is because, in the massless theory, the light cone expansion is given by: 

ZERO < Q 4 > <Q 6 > nrm 

0NP = o h C 4 4 ^ C 6 fi 1 ( 107 ) 

m~ ntf. 

In fact there are no dim-2 Lorentz and gauge invariant operators. For example, is not gauge invari- 
ant. In the massive theory, the ZERO is replaced by light quark mass-squared m 2 . This is still negligible 



if m is taken as a lagrangian mass of a few MeV. But would not at all be negligible, actually would very 
much affect the result, if it is taken as a constituent mass of order m ~ Aqcd- Most people believe the 
optimistic version. I am not convinced that the gap is not filled up by ambiguities of 0(AQ CD /m 2 ) from 
Spert. In any case, one can discuss the error, but it is true and remarkable, that the central value from r 
decay, obtained at very small Q 2 , is in perfect agreement with all other precise determinations of a s at 
more typical LEP values of Q 2 . 



6.2 a s from Deep Inelastic Scattering 

QCD predicts the Q 2 dependence of F(x, Q 2 ) at each fixed x, not the x shape. But the Q 2 dependence is 
related to the x shape by the QCD evolution equations. For each x-bin the data allow to extract the slope 
of an approximately stright line in dlogF(x, Q 2 )/dlogQ 2 : the log slope. The Q 2 span and the precision 
of the data are not much sensitive to the curvature, for most x values. A single value of Aqcd must be 
fitted to reproduce the collection of the log slopes. For the determination of a s the scaling violations 
of non-singlet structure functions would be ideal, because of the minimal impact of the choice of input 
partem densities. We can write the non-singlet evolution equations in the form: 

d i t?i *\ a ^ f 1 dyF(y,t) x 

—logF(x,t) = — — / — — — -P qq {-,a s (t ) (108) 
at 2tt J x y F(x,t) y 



where P qq is the NLO splitting function. It is clear from this form that, for example, the normalisation 
of the input density drops away, and the dependence on the input is reduced to a minimum (also there 
is a single density, while in general there are quark and gluon densities). Unfortunately the data on 
non-singlet structure functions are not very accurate. If we take the difference of data on protons and 
neutrons, F p — F n , experimental errors add up in the difference and finally are large. The F^n data are 
directly non-singlet but are not very precise. A determination of a s from the CCFR data on F^ v n has 



led to [12] 



a s {m z ) = 0.118 ±0.006 (109) 

When one measures a s from scaling violations on F2 from e or \i beams, the data are abundant, the errors 
small but there is an increased dependence on input parton densities and especially a strong correlation 
between the result on a s and the input on the gluon density. There are two most complete and accurate 



derivations of a s from scaling violations in F2. In the first analysis, by Santiago and Yndurain Q13Q, the 
data on protons from SLAC, BCDMS, E665 and HERA are used with NLO kernels plus the NNLO first 
few moments. The analysis is based on an original method that uses projections on a specially selected 
basis of orthogonal polynomials. The quoted result is given by: 

a s {m z ) = 0.117 ±0.003 (110) 



A different analysis by Alekhin [ 14] of a similar collection of proton data from SLAC, BCDMS, NMC 



and HERA with NLO kernels and a more conventional method leads to 

a s {m z ) = 0.116 ±0.003 (111) 
In both analysis the dominant error is theoretical and could be somewhat larger than quoted. 



If we compare these results on a s from DIS with the findings at the Z, given by eq.(103), we see 
that the agreement looks perfect. But, in my opinion, the situation of a s from DIS is not yet completely 
satisfactory (while it is so for a s from e + e~). The data have shown large fluctuations in the recent past. 
For example, the original result obtained by combining BCDMS and SLAC data was [|il|]: 

a s {m z ) = 0.113 ±0.005 (112) 

This low value appeared to be confirmed by the CCFR result from F2 and F3 data combined: 



a s (m z ) = 0.111 ±0.005 



(113) 



But later the same data were corrected by a new energy calibration performed by the collaboration and 
the result become: 

a s {m z ) = 0.119 ±0.005 (114) 

So the experimental errors are perhaps larger than quoted and there is a problem for matching the system- 
atics of different experiments. On the theory side the analysis methods are perhaps still not completely 
optimised. 

Summarising: there is very good agreement among many different measurements of a s (see 
Fig. 20, This is a very convincing, quantitative test of QCD. The average value quoted by PDG 
2000 is 

a s (m z ) = 0.118 ±0.002 (115) 
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Fig. 20: The PDG (Ref. M\) compilation of the existing precise measurements of a s (M2 



The corresponding value of A (for rif = 5) is: 

A 5 = 209 ±24MeV (116) 

A is the scale of mass that finally appears in massless QCD. It is the scale where a s (A) is of order 1. 
Hadron masses are determined by A. Actually the p mass or the nucleon mass receive little contribution 
from the quark masses (the case of pseudoscalar mesons is special, as they are the pseudo Goldstone 
bosons of broken chiral invariance). Hadron masses would be almost the same in massless QCD. 

7. Conclusion 

We have seen that perturbative QCD based on asymptotic freedom offers a rich variety of tests and we 
have described some examples in detail. QCD tests are not as precise as for the electroweak sector. But 
the number and diversity of such tests has established a very firm experimental foundation for QCD as a 
theory of strong interactions. 

The field is still very much in movement. There are areas of continuous development: I list here 
some of these domains where work is in progress and which would deserve an expanded discussion. 

Higher order calculations: important examples are the effort to complete the 3-loop analytic deter- 
mination of the splitting functions, and the calculation of 4-jet distributions at NLO in e + e~ annihilation. 



Resummations in problems with 2 (or more) different large scales: if y=scalel/scale2 one needs to 
resum (a s log a y) n , with a = 1, 2, at all orders in n. We have seen an example on the W, Z pr distribution 
in pp or pp collisions for Aqcd « Pt « mw where all terms of order (a s (pj) logm^/pj,) 71 must 
be taken into account. Other examples involve structure functions at small x ((a s (Q 2 ) log l/x) n ) or at 
x ~ 1 ((a s (Q 2 ) log 1/(1 — x)) n ), or thrust distributions near T ~ 1 etc. 



Renormalons and power suppressed corrections (for an introduction see, for example, Ref Q16|]).. 
The QED and QCD perturbative series, after renormalisation, have all their coefficients finite, but it is 
well known that the expansion does not converge. Actually the perturbative series is not even Borel 
summable. After Borel resummation for a given process one is left with a result which is ambiguous 
by terms typically down by — exp (—n/ba), with n an integer and b the first (3 function coefficient. In 
QED these corrective terms are extremely small and not very important in practice. On the contrary 
in QCD a = a s (Q 2 ) ~ l/(MogQ 2 /A 2 ) and the ambiguous terms are of order (1/Q 2 ) n , that is are 
power suppressed. The problem arises of the precise relation between the ambiguities of the perturbative 
expansion and the higher twist corrections, which is a very interesting field of research. 

New areas explored by experiment. Particularly interesting physics cases are present in polarized 
deep inelastic scattering, or in the phenomenology of structure functions at small x as measured at HERA. 

An important domain of activity is the QCD event simulation for the preparation of LHC experi- 
ment which poses highly non trivial theoretical and technical problems. 

8. APPENDIX: The Formalism of Gauge Theories 

We summarize here the definition and the structure of a gauge Yang-Mills theory. 

Consider a Lagrangian density £[0,9^] which is invariant under a D dimensional continuous 
group of transformations: 

<// = U{e A )<j) (A = 1,2,...,D) . (117) 
For 6 A infinitesimal, U(6 A ) = l + igT,A A T A , where T A are the generators of the group V of trans- 



formations (117) in the (in general reducible) representation of the fields 4>. Here we restrict ourselves 
to the case of internal symmetries, so that T A are matrices that are independent of the space-time co- 
ordinates. The generators T A are normalized in such a way that for the lowest dimensional non-trivial 
representation of the group Y (we use t A to denote the generators in this particular representation) we 
have 

tr(t A t B ) = U AB . (118) 
The generators satisfy the commutation relations 

[T A ,T B ]=iC A B C T c . (119) 
In the following, for each quantity V A we define 

V = T A V A . (120) 

A 

If we now make the parameters 9 A depend on the space-time coordinates 9 A = 9 A (x fl ), C[4>, <9 M </>] is in 
general no longer invariant under the gauge transformations U[6 A (x /1 )], because of the derivative terms. 
Gauge invariance is recovered if the ordinary derivative is replaced by the covariant derivative: 

D, x = d lx + ig\ lx , (121) 

where V A are a set of D gauge fields (in one-to-one correspondence with the group generators) with the 
transformation law 

V; = UV^U- 1 - {l/ig^U'W- 1 . (122) 



For constant 9 , V reduces to a tensor of the adjoint (or regular) representation of the group: 

V; = UV^U- 1 ~ + ig[9, V M ] , (123) 

which implies that 

= v£ - g c ABC e A \^ 

where repeated indices are summed up. 



K° = - 9Cabc0 A V* , (124) 



As a consequence of Eqs. ( |121| ) and (|122j), -D^0 has the same transformation properties as eft: 

(D^)' = U(D^) . (125) 

Thus C[(j), Dn4>] is indeed invariant under gauge transformations. In order to construct a gauge- 
invariant kinetic energy term for the gauge fields V A , we consider 

[Dp, D v ]<j> = igid^Vv - d v \ „ + ig[Vp, V v ]}0 = igF^ , (126) 

which is equivalent to 

< = - ^ - gCABcV*vc . (127) 



From Eqs. ( 1 17 ), ( 125 ) and (|126|) it follows that the transformation properties of F A are those of a tensor 



/if 1 

of the adjoint representation 

K» = tflV^ -1 • (128) 

The complete Yang-Mills Lagrangian, which is invariant under gauge transformations, can be written in 
the form 

£ ym = ~\ E F^F AliU + * D^] . (129) 

For an Abelian theory, as for example QED, the gauge transformation reduces to U[6(x)] = 
exp[ieQ9(x)], where Q is the charge generator. The associated gauge field (the photon), according to 
Eq. ( |122D , transforms as 

V^ = Vp- dp9{x) . (130) 

In this case, the F^ tensor is linear in the gauge field so that in the absence of matter fields the theory 

F 

«n that thp thpnrv i« nnn-trivinl pvpn in thp nh 



is free. On the other hand, in the non-Abelian case the F^ v tensor contains both linear and quadratic 
terms in V/f , so that the theory is non-trivial even in the absence of matter fields. 
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